Two-dimensional electronic transport on the surface of 3D topological insulators 
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We present a theoretical approach to describe the 2D transport properties ol the surfaces of three 
dimensional topological insulators (3DTIs) including disorder and phonon scattering effects. The 
method that we present is able to take into account the effects of the strong disorder-induced carrier 
density inhomogeneities that characterize the ground state of the surface of 3DTIs, especially at low 
doping, as recently shown experimentally. Due to the inhomogeneous nature of the carrier density 
landscape, standard theoretical techniques based on ensemble averaging over disorder assuming a 
spatially uniform average carrier density are inadequate. Moreover the presence of strong spatial 
potential and density fluctuations greatly enhance the effect of thermally activated processes on 
the transport properties. The theory presented is able to take into account all the effects due to 
the disorder-induced inhomogeneities, momentum scattering by disorder, and the effect of electron- 
phonon scattering processes. As a result the developed theory is able to accurately describe the 
transport properties of the surfaces of 3DTIs both at zero and finite temperature. 
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I. INTRODUCTION 

In strong three dimensional topological insulators 
(3DTIs) the nontrivial topology of the bulk energy 
bands 1-3 enforces the essential existence of 2D metallic 
surface states that can be well described at low energies 
as massless Dirac fermions. The valence and conduc- 
tion band of the surface states touch at isolated points, 
the Dirac points, whereas the bulk states are gapped. 
Angle resolved photo-emission spectroscopy experiments 
have confirmed the existence of the Dirac like surface 
states in Bii^SbA Bi 2 Se^, Bi 2 Te^ and Sb 2 Te^. 
Electronically, the surface of a 3DTI is very analogous 
to graphen o 9 ' 10 in which the fermionic excitations are 
also well described, at low energies, as massless Dirac 
fermions. There are two fundamental differences between 
graphene and the surfaces of 3DTIs: (i) In graphene the 
chiral nature is due to the locking of the momentum di- 
rection with the electron pseudospin, associated with the 
sublattice degree of freedom, instead of the real spin as 
in the surface of strong 3DTIs; (ii) In strong 3DTIs the 
number of Dirac points is odd whereas in graphene is 
even. These differences make graphene and the surfaces 
of strong 3DTIs fundamentally different. However, the 
fact that graphene and the surfaces of strong 3DTIs are 
both two-dimensional electronic systems and have a very 
similar band structure suggests that these two systems 
might have similar charge-transport properties. As we 
show in this work this is only partially correct. 

Two important aspects of TI surface transport need 
to be mentioned (in the context of our comprehensive 
theoretical work to be presented in this paper) so as to 
avoid any confusion about our goal and scope. First, 2D 
TI transport occurs on the surface of 3D TI materials, 
and theoretically the bulk 3D states should be insulating 
with no metallic contribution to the conductivity. Exper- 
imentally, however, this situation of 2D metallic trans- 
port on a bulk 3D insulator has not yet been achieved 



since the bulk, instead of being a band insulator, seems 
to have a lot of free carriers which contribute (indeed, 
often dominate) the measured conductivity^. We com- 
pletely ignore the complications of the bulk conduction 
and the bulk carriers in our theory concentrating entirely 
on purely the 2D surface conductivity (as a function of 
density and temperature) assuming the bulk to be an in- 
sulator as the theory implies it should be. Recent exper- 
imental advances in materials preparation and thin film 
device fabrication have made it possible to see the ex- 
pected pure 2D surface conduction with little contamina- 
tion from the bulk states, and thus, our work is relevant 
to an increasing body of recent data^— . In any case, 
the subject is interesting only because of the 2D metallic 
surface topological states, and therefore, we focus en- 
tirely on this issue. Second, the actual energy dispersion 
of the 3DTI surface states appears to follow the linear 
Dirac-like spectrum only at rather low energies, partic- 
ularly for the hole states (i.e. the valence band) with 
strong nonlinearity becoming apparent at higher ener- 
gies. This nonlinearity (e.g. a quadratic correction to the 
linear dispersion) is nonuniversal and strongly materials- 
dependent whereas the theoretical Dirac-behavior at low 
energies is universal (with different Fermi velocities char- 
acterizing different TI materials). Because the focus of 
our work is on the universal behavior arising from den- 
sity inhomogeneity (which is much more important at 
low Fermi energies, i.e. low 2D densities), which should 
not depend much on the details of the band structure at 
higher energies, we ignore the higher energy parabolicity 
of the 2D surface bands. We emphasize that these two 
approximations, neglect of bulk conduction and surface 
band dispersion correction at higher energy, imply that 
our theory should not be construed as a quantitative the- 
ory for any particular TI material, but as a qualitative 
guide for the universal features of 3DTI surface transport 
behavior. Detailed quantitative comparison between the- 
ory and experiments is further complicated by our lack of 
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knowledge of the precise parameters for TI systems any 
way (e.g. the Fermi velocity, the phonon parameters, the 
nature of disorder, etc.), and therefore, our theory pro- 
vides the zeroth order qualitative theory for 3DTI surface 
transport which should apply to all TI materials. 

The surfaces of the newly discovered strong 3DTIs are 
of great fundamental interest*^— and in addition have 
the potential to be used in disruptive novel technologies 
such as topological quantum information processing^ - - — . 
To be able to use the surfaces of 3DTIs to study novel fun- 
damental phenomena and in novel technologies it is es- 
sential to understand their electron transport properties 
and in particular to understand the main factors limit- 
ing their electron mobility. Recently, experiments on thin 
films of strong 3DTIsi2~— , by enhancing the surface-to- 
volume ratio^ 3 -, have been able to greatly reduce the bulk 
contribution transport and obtain the intrinsic 2D trans- 
port properties of the surfaces of strong 3DTIs. Previous 
theoretical works on the transport in the surface of strong 
3DTI s 24 ' 25 have used simplified models and have mostly 
focused on the zero temperature limit. 

In this work we present a comprehensive transport the- 
ory for the surfaces of 3DTIs valid both at zero tempera- 
ture and at finite temperatures. As in graphen o 26 ' 2 , one 
of the main difficulties in developing a transport theory 
for the surfaces of 3DTIs is the presence, especially at 
low doping, of strong carrier density inhomogeneities^ 
induced by disorder. The presence of strong spatial fluc- 
tuations makes the standard theoretical approaches, that 
rely on the homogeneous nature of the ground state, inad- 
equate. Moreover, because of the inhomogeneities, at fi- 
nite temperatures the contribution of thermally activated 
carriers to transport can be very important and difficult 
to quantify. Thus, disorder-induced density inhomogene- 
ity by itself could introduce considerable insulating-like 
activated transport behavior in the nominally metallic 
3DTI surface conduction. The transport theory that we 
present takes into account both the effects of quenched 
disorder and electron-phonon scattering processes. In 
particular our theory is able to take into account the ef- 
fects of the strong disorder-induced inhomogeneities both 
at zero and finite temperature. 

To characterize the inhomogeneous ground state we 
generalize to the case of 3DTIs' surfaces the Thomas- 
Fermi-Dirac- Theory (TFDT) first developed to study 
graphen o 10 ' 29 . Combining the TFDT results and the 
Boltzmann theory we develop and validate the effective 
medium theory (EMT) to obtain the transport proper- 
ties at zero temperature. We then develop an effective 
2-fluid transport theory that we validate at zero temper- 
ature by comparing its results to the ones obtained us- 
ing the EMT. The great advantage of the 2-fluid theory 
is that it allows to readily obtain the transport proper- 
ties at finite temperature including all the temperature 
dependent effects: electron-phonon scattering, thermal 
activation, changes with temperature of the screening 
properties, and thermal broadening of the Fermi sur- 
face. Given the experimental evidenc o 12 ' 30 that in cur- 



rent 3DTIs charged impurities are the dominant source of 
disorder we have applied the theory to the case in which 
the quenched disorder is due to random charges placed 
in the vicinity of the surface of the 3DTI. We provide 
both a detailed analytical theory and comprehensive nu- 
merical results for the density and temperature depen- 
dent 3DTI surface transport properties in the presence 
of density inhomogeneity, scattering by random charged 
impurities, and phonon scattering. 

In section [II] we present the theoretical approach that 
we have developed to describe the transport on the sur- 
faces of 3DTIs taking into account both quenched disor- 
der and electron-phonon scattering processes, in particu- 
lar when the quenched disorder is due to charged impuri- 
ties. In section Ull Al we present our results for the char- 
acterization of the disorder-induced carrier density inho- 
mogeneities obtained using the TFDT. In section llll Bl wc 
present our results for the conductivity at zero tempera- 
ture and finally in section llll Cl we present our results for 
the conductivity at finite temperature. Section HVl briefly 
summarizes our findings and the differences between the 
2D transport properties of the surfaces of 3DTIs and sin- 
gle layer graphene. 



II. THEORETICAL APPROACH 

To study the electronic transport on the surface of 
strong 3DTIs we use the Boltzmann theory. From the 
Boltzmann theory, within the "relaxation time approx- 
imation" , the electronic conductivity a is given by the 
following equation: 

a(n,T) = deD(e)v F r tot (e,T) (-^) (1) 

where e is the electron charge, D(e) the density of states 
at energy e, vp the Fermi velocity, /(e) = 1/(1 + e^^ e ~^) 
the Fermi-Dirac distribution, and r tot the total trans- 
port mean free time due to the electron scattering off 
quenched disorder and phonons. Assuming independent 
scattering from disorder and phonons (we mention that 
this is not equivalent to assuming the Matthiessen's rule 
which assumes that the resistivity due to independent 
mechanisms can be added and is invalid for our system) 
we have 



7~tot ^~dis Tpfa 

where Tdis is the transport mean free time due to elec- 
trons scattering off quenched disorder and r p h is the 
transport mean free time due to electron-phonon scat- 
tering processes. One thing that we must emphasize is 
that Eq. (Q} is valid as long as the system is homogeneous 
(i.e. spatial density fluctuations effects are small enough 
so that the average density n is a meaningful quantity, 
an approximation which would break down for low n). 
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The energy dependent scattering time Tdj S (e) due to 
quenched disorder is given by 



where 



7~dis (^pk) 



= 2irrid, 



d 2 k' 



|(^k,pk<)| 2 2(0kk<) 



x [1 - cos 0kk'] £(e P k' - e P k) (3) 



where e pk = phvp\k-\ is the energy of a quasiparticle with 
p = ±1 and momentum k, n^is is the 2D density of impu- 
rities, (Vpk.pk') is the matrix element of the scattering po- 
tential, g(0kk>) = [1 + cos6*kk'] /2 is the TI chiral matrix 
element factor arising from the wavefunction overlap be- 
tween states with momentum k and momentum k' with 
#kk' the angle between k and k'. In Eq. ([3]), to minimize 
the number of parameters entering the theory we have 
assumed that the impurities are randomly distributed in 
a 2D plane located at an effective distance d from the sur- 
face of the 3DTI. It is straightforward to include in the 
theory a more complex three-dimensional distribution of 
quenched impurities, but given the lack of experimental 
information about the distribution of unintentional and 
unknown quenched impurity disorder in the system, it 
is theoretically more meaningful to use a minimal model 
with just two unknown parameters Udts and d, which can 
simulate essentially any realistic disorder distribution in 
an approximate manner- we note that d = implies that 
the charged impurities are simply located on the surface 
of the 3DTI. 

There is considerable evidence 1 ^ that in 3DTIs ran- 
dom unintentional charged impurities are the dominant 
source of disorder scattering. We therefore assume that 
the quenched disordered potential Vd is due to charged 
impurities. For charged impurities, taking into account 
screening by the 2D surface carriers themselves in addi- 
tion to the screening by the background lattice, we have 
ndis\{V p u, P k>)\ 2 = n imp \vi{q)/e{q,T)\ 2 where v l {q) = 
2ne 2 /(nq) is the Fourier transform of the 2D Coulomb 
potential in a medium with an effective background static 
lattice dielectric constant k, and e(q, T) is the 2D static 
RPA dielectric function at finite temperature 3 ^. The re- 
ported values of k for Bi2Se3 range from 30, Ref. [28|, to 
« 55, Ref. GH. 

In Bi2Se3 the lowest optical phonon energy has been 
measured to be 8.94 meV~ 100 K 32 i 33 and therefore the 
optical phonons provide a substantial source of scatter- 
ing only at high temperature (> 250 K)^. Because we 
are interested primarily only in the transport properties 
at temperatures below 250 K in the remainder we ne- 
glect the contribution to the resistivity due to optical 
phonons and consider only the contribution due to acous- 
tic phonons. 

Following references [351 - 4371 we have that, considering 
only longitudinal acoustic phonons, the scattering time 
r p h is given by 



T p h(e) 



B 1 



cos0 kk /)W kk . 



W(fO 
W(£)' 



9-ir 

ww = T i^q)! 2 x W( £ - £ '+ ^) 

q 

+ (N q + l)6(e-s'-^)], (5) 

is the transition probability from the state with momen- 
tum k to the state with momentum k'. In {2} q = k — k', 
C(q) is the matrix element for scattering by acoustic 
phonon, cj q = t>;q is the acoustic phonon frequency with 
vi the phonon velocity, and N q — l/(exp(/3w q ) — 1) is the 
phonon occupation number. The matrix element C(q) 
for the deformation potential electron-phonon coupling 
is given by 



|C(q)| 2 



D 2 hq 
2Ap m vi 



1 



\2k) 



(6) 



(4) 



where D is the deformation potential coupling constant, 
A is the area of the sample, and p m is the 2D mass density 
of one quintuple layer (around 1 nm thick) of Bi2Se3, 
given that the length scale over which the 2D surface 
states decay into the bulk is approximately 1 nmas. 

Using Eqs. ((H)-® we can calculate the conductiv- 
ity taking into account scattering events due to both 
quenched disorder and phonons as long as the system 
is homogeneous. However, especially close to the Dirac 
point, the random charged impurity induced disorder po- 
tential causes the carrier density landscape to become 
strongly inhomogeneous, a fact that has been observed 
experimentally in TIs 2 ^ and previously in graphen o 26 ' 27 . 
To develop a theory in the presence of strong inhomo- 
geneities it is first necessary to characterize them. To 
do this we use the Thomas-Fermi-Dirac- Theory (TFDT) 
first introduced in Ref. 29; In the TFDT, similarly to 
the Density Functional Theory (DFT), the energy of the 
system is given by a functional of the density profile n(r). 
The great advantage of a functional formalism is that it 
is not perturbative with respect to the spatial fluctua- 
tions of the carrier density and therefore can take into 
account nonlinear screening effects that dominate close 
to the Dirac point. TFDT is just well-suited to describe 
the situation with large disorder-induced spatial density 
inhomogeneity as in the low carrier density case whereas 
in the high-density situation, it simply gives the homoge- 
neous density result with small fluctuations around the 
average density. In the TFDT, contrary to DFT, also 
the kinetic energy term is replaced by a density func- 
tional. This simplification makes the TFDT very effi- 
cient computationally and therefore able to obtain disor- 
der averaged quantities, a task that cannot be accom- 
plished using DFT. The simplification also makes the 
TFDT in general less accurate than DFT—, however as 
long as the characteristic length-scale over which the den- 
sity varies is larger than the local Fermi wavelength X F , 
i.e. \Vnfn\- 1 » \ F 2Z£2£L the TFDT returns reliable 
result s 29 ' 42 ' 43 . Our results show that as in graphene 29 
the condition |Vn/n| _1 ^> \p is satisfied for the sur- 
face of 3DTIs in typical experimental conditions. Close 
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to the charge neutrality point (CNP) the density inside 
the electron-hole puddles is always different from zero 
(so that Xf is always finite) and of the order of n rms . 
Our results show that n rms ~ rii mp and therefore the 
TFDT is also valid at the CNP as long as rii mp is not 
too small. The great advantage of TFDT over DFT (to 
which TFDT is an approximation, as it uses the non- 
interacting kinetic energy functional) is that its relative 
numerical and computational ease enables one to use it 
for the calculation of transport properties using the com- 
puted ground state inhomogeneous spatial density pro- 
file, which would be completely computationally impos- 
sible for DFT to do. 

Using the TFDT we can characterize completely the 
carrier density profile in the presence of a disorder poten- 
tial. We can obtain the typical length scale L d and root 
mean square fluctuation n rms of the disorder- induced car- 
rier density inhomogeneities. Using the Boltzmann the- 
ory we obtain the relation between the mean free path 
and the doping, l(n) — vpTt t(n), valid in the homoge- 
neous limit. In the limit in which £{n rms ) -C Lp the 
number of scattering events inside a single homogeneous 
region, puddle, of the inhomogeneous landscape is large 
enough that the Boltzmann theory is valid locally. In ad- 
dition, due to the Klein tunneling, as in graphene^~— , 
the resistance due to the boundaries between the puddles 
can be neglected in comparison to the resistance arising 
from scattering events inside the puddles^. Under these 
conditions, due to the random distribution of the pud- 
dles, 2D transport on the surface of a 3DTI can be de- 
scribed by the effective medium theory (EMT) 10 i 40 i 51 ~.ii. 
In the EMT, which is extensively used in science and en- 
gineering to quantitatively describe properties of highly 
inhomogeneous systems, the conductivity of the inhomo- 
geneous system is obtained as the conductivity ctemt 
of an equivalent homogeneous effective medium by aver- 
aging over disorder realizations the local values a(n(r)) 
given by the Boltzmann theory. The resulting implicit 
equation for oemt is: 



dn 



cr(n) - uemt 



cm 



oemt 



P[n] = 



(7) 



where P[n] is the disorder-averaged carrier density prob- 
ability distribution that we obtain using the TFDT. A 
solution of the implicit EMT integral equation defined by 
Eq.[7]providcs the effective conductivity of the inhomoge- 
neous system. In graphene the TFDT+EMT method has 
been shown to give results in remarkable agreement with 
experiments^ - even in the highly inhomogeneous situa- 
tion very close to the Dirac point and with full quantum 
transport analysis 50,55 . 

A simplified approach that allows to make further an- 
alytical progress and obtain results in qualitative agree- 
ment with the numerical TFDT-EMT approach is the "2- 
fluid" model 5 ^. In this approach the inhomogeneous state 
characterized by the presence of electron-hole puddles is 
approximated as a system comprised of the "electron- 
fluid", formed by the electrons, and the "hole-fluid" 



formed by the holes, with conductivity a e and ah re- 
spectively. Let p be the fraction of the system occupied 
by the electron gas, and consequently (1 — p) the fraction 
occupied by the hole gas. Adapting Eq. ([7]) to the sim- 
ple case of only two components for the inhomogeneous 
system we obtain the effective conductivity 5 ^—: 



1. 



(cr e - a h ) + 



<J h f 



{2p-lf 



(8) 

One advantage of the 2-fluid model is that it allows to 
easily take into account the effect of activation processes 
that at finite temperature, especially close to the Dirac 
point, give a substantial contribution to the conductivity 
and in particular qualitatively modify its temperature de- 
pendence. Carrier activation becomes operational when 
local potential fluctuations due to weak screening at low 
density lead to carrier confinement or localization in pud- 
dles, and global transport involves thermal activation of 
carriers over the local potential hills and barriers. The 
activation process, which obviously becomes more im- 
portant as the inhomogeneity becomes more important 
at lower carrier density, cannot be captured by the sim- 
ple Boltzmann theory of Eq. (JXJ) or, for that matter, by 
any ensemble-averaged transport theory. To take into ac- 
count the presence of activation processes a e and <7h can 
be written as a sum of two terms 5 ^: 

a e = + 4 a) 

*h = 4" a) + 4 a) 

where ai na ^ , <J^°^ are the disorder averaged conductivi- 
ties obtained using the Boltzmann theory and Oe , 
the contribution to the conductivity from activation pro- 
cesses. 

<7e j ( cr i" a '') can be obtained by multiplying 
Eq. (HJ) by the ratio n e /n (n/j/no) 56 - Here, n e — 
JZ D e (e)f(e)de (n h = JZo^i 1 ~ denotes 
the effective electron (hole) density of inhomogeneous 
systems, while no — f Q °° 2 ^(hvl)' 1 /( e )^ £ * s the electron 
density of homogeneous systems. The density of states 
D e (e) (Dh(e)) after disorder averaging are given by££>: 



D e (e) 
D h {e) 



g s gv(e-V) 



, 2ir(hv F ) 2 
9s9v(V~e) 
27t(/juf) 



P(V)dV 
2 P(V)dV 



where P(V) is the probability distribution of the screened 
disorder potential. We want to mention that we use the 
density of states for homogeneous systems in Eq. ((T|), 
i.e., D(e) = 2 ^(hvl) ,i ' ^° av °id double counting since the 
density inhomogeneity effects have already been consid- 
ered through the variation of effective carrier density. At 
higher doping, both n e /riQ and the fraction of area oc- 
cupied by electrons, denoting as p = J_ F P(V)dV, ap- 
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proach unity while Uh/na approaches zero. Note, how- 
ever, that the ratio n e /n (rih/no) being temperature- 
dependent is generally not equal to p ((1 — p)). Given 
that the 2-fluid model is an effective model the use of the 
exact P(V) does not guarantee an increase of its accu- 
racy. It is more sensible to simply assume P(V) to have 
an effective Gaussian profile 



P(V) 



1 



V2 



ITS' 



: exp(-V 2 /2s 2 ). 



(9) 



with effective variance s 2 . The TFDT results, presented 
in Sec. IIII1 show that the root mean square s of the 
screened disorder potential depends weakly on the dop- 
ing or carrier density. In Eq. [5] we can then neglect the 
dependence of s on the average doping and use the value 
obtained for the Dirac point. We consider two ways to 
estimate the value of the effective s at the Dirac point 
that enters the Gaussian approximation for P(V): (i) the 
self- consistent approximation introduced in Ref. l59t (ii) 
the quasi- TFDT approximation in which s is fixed using 
the relation between rt rms and s in the Thomas-Fermi 
approximation with the value of n rms obtained from the 
full TFDT calculation. 

In general P(V) cannot be obtained analytically 
but it is possible to obtain explicit expression for its 
moments 60 . For the second moment ((SV) 2 ) taking 
into account screening effects, within the random-phase- 
approximation (RPA) for the surface states of a 3DTI 
with total degeneracy g = g s g v = 1 we have 



,60. 



((SV) 2 ) = 2nn lmp 
where 

e -4a 

C (r s ,a) = -l + i— 



C (r s ,a=k F d) (10) 



64^ [4a] 



(11) 

E l [Aa + 2r s a]), 



+ e 2r ° a (l + 2r s a)(E 1 [2r s a] 

r s = e 2 /(hv F n), and Ei[z] = t~ 1 e~ t dt is the expo- 
nential integral function. Assuming P(V) to be a Gaus- 
sian with variance s 2 we have s 2 = ((SV) 2 ). The diffi- 
culty arises from the fact that the function Co depends 
on the density via the Fermi wavevector k F that, in the 
presence of disorder, at the Dirac point cannot, because 
of fluctuations, be taken to be simply zero. In the self- 
consistent approximation one assumes that at the Dirac 
point the system can be approximated by a homogeneous 
system having an effective carrier density h = (g/4ir)kp 
such that Ep = s 2 , i.e. H 2 vpkp = s 2 . Using this rela- 
tion and Eq. (JTUJ) we obtain the following self-consistent 
equation for s: 



2nnr 



C (r s ,a = sd/(hv F )). 



(12) 



By solving Eq. (|12l) we obtain the value of s within the 
self-consistent approximation. In the self- consistent +2- 
fluid model the variance of the effective Gaussian proba- 
bility distribution P(V) is obtained using Eq. (fT2")l . 



By minimizing the Thomas-Fermi energy functional, 
at the Dirac point, for g = 1, we obtain the following 
equation 

2hv F sgn(n) ( ^ n ^ \ ' _ y(r) = (13) 
V 9s9v J 

where the first term is due to the kinetic energy and 
V(r) is the local value of the screened disorder potential. 
By disorder averaging (fT31 and assuming P(V) to be a 
Gaussian with variance s 2 we obtain 



9v \/3 



2-k(Hvf) 2 4 



(14) 



In the quasi- TFDT approximation the value of s used in 
the 2-fluid model is set using Eq. (jT4")) and the value of 
n rms obtained from the TFDT at the Dirac point. 

Locally, the activated conductivities ai a \ a^ are 
given by: 

a^ (V(v)) = *M exp[/3(£ F - V(v))} (15) 
al a) (V(v)) = exp[/J(V(r) - E F )]. (16) 

where V(r) is the local value of the screened disorder 
potential. By disorder averaging these expressions and 
summing the contribution of the non-activated conduc- 
tivities we finally find 56 : 



Oh 



i/ 

P J- 


■ E f 

-oo 


(a (na) . 




■>4 


(na) 

■ i « 




2p 


i 




p OO 

/ 


(i- 


P) 


4 na 


' + 


(na) 




2(1-3 



E F 



+ 



Ps\ 



V2s V2J 



(.17) 



(na) 
h 



r^)P(V)dV 



^+^erfc 



/3s 
T2 



18) 



where the fraction of the system occupied by electrons is 
given by p = j_ F P(V)dV with Fermi energy E F . These 
are the expressions of the electron and hole conductivi- 
ties that enter the 2-fluid model whose total conductivity 
at is then obtained by using Eq. (|8|). We notice that in 
the expressions (fl7|) . (fT8|) several temperature effects are 
taken into account: (i) the effect due to clcctron-phonon 
scattering processes that affect the total scattering time 
and therefore the value of ai^ and a h na ^; (ii) the tem- 
perature dependence of the dielectric functions that en- 
ters in the calculation of Tdis and that affect the value 
of ai na ^ and a^°^] (iii) the thermal broadening of the 
Fermi surface; (iv) the presence of thermal activation. 
The ability to capture all these effects makes the 2-fluid 
model a very useful tool to study the transport properties 
of disordered TIs. 

It is useful to explicitly write down the expression for 
the conductivity acNP at the charged neutrality point 
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obtained using the 2-fluid model. At the CNP we have 
that 50% of the sample will be covered by electron pud- 
dles and 50% by hole puddles so that p = 1 — p = 0.5. In 
this case at = a e = <7h- At zero temperature, there are 
no thermal activation effects and we find£&: 



&CNP — cr, 



(na) 
CNP 



1 



9s9v 



8G[r s /2] h 2ir(hv F ) 2 An mp 



where 
G[x 



3ttx 2 x(3x 2 — 2) arccos[l/a;] 
Using Eq. (fT4"]) we can rewrite (fK)f in terms of n rm 



,, = -+3x 
2 4 



(na) 

T CNP 



1 



8G[r s /2] h 



(19) 



(20) 



(21) 



In the following sections we present our results and 
in particular, when possible, a comparison between the 
three approaches introduced: TFDT+EMT approach, 
quasi-TFDT+2-fluid model, and self-consistent+2-fluid 
model. 



III. RESULTS 

A. TFDT results 

In this section we present the results for the carrier 
density distribution. These results will then be used to 
calculate the conductivity within the EMT and the 2- 
fluid models. 

In Fig. [TJ we show the carrier density profile calculated 
using the TFDT for a single disorder realization with 
charged impurity density n imp — 7.5 x 10 13 cm -2 and two 
different values for the impurity distances, d = 0.1 nm, 
and d = 0.2 nm in panel (a) and (b) respectively. Un- 
less otherwise specified, we use the background dielec- 
tric constant n = 50 for Bi2Se3 and the Fermi velocity 
6.4 x 10 5 m/s!2i£!. 
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FIG. 1. (Color online). Color plots of carrier density distri- 



bution n(r) at the Dirac point for m 



7.5 x 10 13 cm - 



and k = 50. The color-scale is in units of 10 cm . (a) The 



impurity distance d = 0.1 nm. 
d = 0.2 nm. 



(b) The impurity distance 



Figure [T] conveys the nature of the carrier density land- 
scape on the surface of disordered 3DTIs as also shown 
recently by direct imaging experiments^,. To be able to 
make a quantitative comparison with the experiments, 
and calculate the conductivity for large samples, using 
the TFDT we calculate the disordered averaged density 
probability distribution P(n) for different parameter val- 
ues like doping, d and k. As in graphene^ we find that 
P(n) obtained using the TFDT is bimodal, especially for 
finite values of the average carrier density, and so it is 
not well fitted by any single curve. To exemplify this 
finding Fig. [5] shows P(n) at the CNP obtained using 
the TFDT and possible fitting curves. It is obvious from 
the figure that a reasonable fit can be obtained only by 
using two different Gaussian curves, one to fit the very 
high peak centered at n = and one to fit the long tails 
of the distribution. 
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FIG. 2. (Color online). Density probability distribution at 
the Dirac point for rii mp — 7.5 x 10 13 cm~ 2 , d = 0.2 nm and 



k = 50. 



The black line shows the TFDT result. 
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The blue 

line is the fit P = be~ ln/a] with a = 0.7 x 10 12 cm -2 and 
b — 0.5309. The red line is the Gaussian fit to the tails of 
P(n), P = be~ x2/a2 with a = 2.5 x 10 12 cm -2 and b = .125. 
The green line is the Gaussian fit to the center of P(n), P = 



be 



-x /a 



with a = 0.25 x 10 12 cm" 2 and 6 = 1.0. 



The knowledge of P(n) allows the calculation of all the 
statistical properties that characterize the strongly inho- 
mogeneous ground state of the surface of a 3DTI in the 
presence of disorder. The quantity that better quanti- 
fies the strength of the carrier density inhomogeneities is 
n rms , shown in Fig. [3] as a function of the average car- 
rier density (n) for different values of n 
larger values of n 
interesting result is that n rms also increases with (n) but 
the dimensionless ratio decreases with increasing average 
density. This is due to the fact that as (n) increases the 
range of values that n can take locally also increases in- 
ducing larger values of n lms , but the dimensionless ratio 
of the density fluctuation to the average density decreases 
with increasing density. 



imp . As expected 
lmp induce larger value of n rms . The 
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FIG. 3. (Color online). Root mean square of the density fluc- 
tuations n rm s as a function of the average density for different 
values of riimp with d = 0.1 nm and re = 50. 



creases. Increasing the average distance of the charged 
impurities from the surface of the 3DTI also strongly re- 
duces the amplitude of the spatial inhomogeneities and 
therefore of n rms and Stfdt as shown in Fig. [5] 
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FIG. 5. (Color online). TFDT results at the CNP for d = 0.1 
nm. (a) n rms as a function of substrate dielectric constant 
re. (b) Root-mean-square of the screened disorder potential 
stfdt as a function of substrate dielectric constant re. The 
solid lines are for rii mp — 4 x 10 13 cm -2 and the dashed-lines 
are for riimp = 11.0 x 10 13 cm -2 . 



Fig. @] (a) shows the TFDT results for the root mean 
square of the screened disorder potential stfdt as a 
function (n) for different values of n lmp . As (n) in- 
creases the screening becomes more effective and there- 
fore stfdt decreases. However, we see that the depen- 
dence of stfdt on (n) is fairly weak, for all values of 
"-imp considered. This fact justifies the assumption in the 
quasi-TFDT model to assume the variance of the effective 
distribution P(V) to be independent of (n). Fig. [4] (b) 
shows the scaling of the ratio stfdt/ y/riimp versus (n). 
We can see that the ratio changes by less than 10% over 
a wide range of experimentally relevant values of the av- 
erage density. This result is one more evidence of the 
weak dependence of stfdt on (n) — to the zeroth or- 
der, the inhomogeneity and fluctuations are determined 
by the impurity distribution rii mp and d (as well as the 
background dielectric constant re). 
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FIG. 4. (Color online), (a) Root mean square of the screened 
disorder potential stfdt as a function of average carrier den- 
sity for different values of rii mp with d = 0.1 nm and re = 50 
using TFDT. (b) stfdt / ' y/n imv (in units of v / 10~ 13 meV-cm) 
as a function of average carrier density. 



The dependence of n rms and Stfdt on re and d is very 
strong as shown by Fig. [5] and Fig. [51 From Fig. [5] we 
see that both n rms and stfdt decrease rapidly as re in- 
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FIG. 6. (Color online), (a) and (b) n rms at the Dirac point 
versus the distance d of the charged impurity to the TI sur- 
face calculated within TFDT for n imp = 4.0 x 10 13 cm -2 . 
The solid, dashed and dot-dashed lines correspond to chemi- 
cal potential Ef = 0, 40 and 100 meV, respectively, (c) and 
(d) The corresponding stfdt results versus d. (a) and (c) 
The effective background dielectric constant re = 5. (b) and 
(d) The effective background dielectric constant re = 50. 

We now present a comparison of the TFDT results with 
the ones obtained using the self-consistent and the quasi- 
TFDT approaches, methods in which P(V) is assumed 
to be Gaussian. Fig. [7] shows the scaling of the screened 
disorder rms obtained using the self-consistent approach 
(s sc ) versus doping and re. These results, analogously to 
the TFDT results, show that s depends very weakly on 
(n) and quite strongly on re. 

Fig. [5] shows the comparison for the value of s at the 
CNP as a function of d obtained using the three different 
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FIG. 7. (Color online). Self-consistent results, (a) s sc / ^/ni mp 
(in units of VlO _13 meV-cm) as a function of carrier density 
for different values of n; mp with k = 50 and d — 0.1 nm. (b) 
Potential fluctuation s sc as a function of substrate dielectric 
constant k. The solid lines are for n,i mp — 4x 10 13 cm -2 and 
the dashed-lines are for n; mj) = 11.0 x 10 13 cm" 2 with d = 0.1 



methods: TFDT, quasi- TFDT and, self-consistent. We 
see that the self-consistent approach in general returns 
values of s larger than the ones obtained using the TFDT. 
We should emphasize that s obtained using the quasi- 
TFDT method is only an effective quantity that is used 
to calculate the transport properties within the 2-fluid 
model. 
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FIG. 8. (Color online). Comparison of three methods for the 
potential fluctuation s at the Dirac point versus the charged 



impurity distance d to the TI surface for m r , 



4.0 x 10 1 



cm" . The solid, dashed and dot-dashed lines correspond 
to the TFDT, quasi- TFDT and the self-consistent methods, 
respectively, (a) The effective background dielectric constant 
K = 5. (b) The effective background dielectric constant « = 
50. 



B. Transport at zero temperature 

Using the TFDT and the EMT we can calculate the 
2D conductivity on the surface of a 3DTI. Fig. [9J shows a 
as a function of doping, i.e. average density, for several 
values of n,i mp and fixed d = 0.1 nm and K = 50. As in 
graphene 1 ^ the TFDT+EMT results recover the behavior 
of <r(n) observed experimentally 1 —: the linear scaling of 
cr(n) at large doping, the finite value (<7 m i n ) of a for n = 0, 
and the crossover regime for intermediate values of n. 
More importantly, the theory returns values of <7 m i n = 



(2 — 4)e 2 /h (depending on the sample properties) that 
agree with the ones observed in experiments 1 ^. 
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FIG. 9. (Color online). TFDT-EMT conductivity as a func- 
tion of average carrier density for different values of rn mp with 
d = 0.1 nm and k = 50. 



Fig. [TU] shows the comparison for the dependence of 
a with respect to n obtained using the TFDT+EMT 
method and the 2-fluid approximations. From Fig. [TUJ 
we see that for nj mp < 4 x 10 13 cm" 2 and values of 
k > 50, the three approaches give very similar results. 
For smaller values of k the three approaches give results 
that are qualitatively similar but that differ quantita- 
tively. The general conclusion is that away from the 
CNP the three approaches agree quantitatively for sam- 
ples with mobility fj, > 1000cm 2 /V • s. For very low mo- 
bility highly-disordered TI systems (with 2D surface mo- 
bility lower than 1000 cm 2 /V-s even at T — ), we expect 
the TFDT to provide the most quantitatively accurate 
results. However, given that our main objective is to de- 
scribe the universal qualities of the transport arising from 
the presence of inhomogeneities, and given the absence 
of accurate knowledge of the impurity distribution, i.e. 
n-imp and d in our model, for the purposes of this work 
the three approaches appear to be equivalent. 

Close to the CNP the three transport approaches give 
results that differ quantitatively as shown in Figs. [TTT 
[T^l The self-consistent+2-fluid model for the parameter 
values relevant for TIs gives values of <r m in smaller than 
the ones obtained using the TFDT+EMT method. In 
TIs the agreement between the two methods close to the 
CNP is worse than in graphene^ 2 - due to the fact that 
in TIs the density of charged impurities is larger than in 
typical graphene samples. It is remarkable to see how, 
in agreement with experiments, the three methods give 
that <7 m in depends very weakly on ni mp and k, for exper- 
imentally relevant parameter values. The weak depen- 
dence on n lmp can be qualitatively understood using the 
2-fluid model result for <7 m i n Eq. (fTTJj) from which we see 
that <7 m i n is proportional to the ratio s 2 /n,i mp , and that 
on the other hand s 2 is proportional to ni mp , Eq. 03 



FIG. 10. (Color online). Calculated conductivity as a func- 
tion of carrier density for different values of the substrate di- 
electric constant with the impurity density rii mp = 4.0 x 10 13 
cm -2 , d = 0.1 nm. The solid, dashed and dot-dashed lines 
are obtained by using TFDT-EMT method, self-consistent 
method, and quasi- TFDT method, respectively. 
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Within the semiclassical approach the very weak depen- 
dence of cr m in on n\ mp and k is due to the fact that an 
increase (decrease) of ni mp (k) increases the strength of 
the disorder potential that causes a decrease of the car- 
riers mean free path £, and an increase of the amplitude 
of the carrier density inhomogeneities(i.e. the density of 
carriers in the electron- hole puddles) . The reduction of £ 
and the increase of n rms have opposite effect on a m - ln and 
they almost cancel ou t 40 ' 63 . 
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FIG. 11. (Color online). The minimum conductivity a m in 
as a function of charged impurity density n; mp with k — 50. 
The symbols diamond, circle and square denote the results 
obtained using the quasi- TFDT, TFDT-EMT and the self- 
consistent method, respectively. The solid and dashed lines 
correspond to the impurity distance d = 1 A and d = 2 A , 
respectively. 



The dependence of a mm with respect to the average 
distance d of the charged impurities from the TI's sur- 
face is appreciably different for the three methods. The 
sclf-consistent+2-fiuid model returns a very weak depen- 
dence of <7 m i n with respect to d. This is due to the weak 
dependence with respect to d of s 2 obtained using the 



FIG. 12. (Color online). The minimum conductivity a min 
as a function background dielectric constant re for rn mv = 
4.0 x 10 13 cm" 2 , d = 0.1 nm. The dashed, solid and 
dot-dashed lines denote the results obtained using the quasi- 
TFDT, TFDT-EMT and the self-consistent method, respec- 
tively. 



self-consistent approximation, Eq. p^|) . On the other 
hand the value of s 2 obtained using the TFDT is quite 
sensitive to the value of d and as a consequence using 
the TFDT+EMT and the quasi-TFDT+2-fluid models 
we find that the dependence of cr m ; n on d is not as weak 
as the one given from the self-consistent +2-fluid model. 
All the three approaches show that cT min decreases as a 
function of d. This is due to the fact that as d increases, 
due to the weakening of the disorder potential, the de- 
crease of n rms is faster than the increase of the mean free 
path. 
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FIG. 13. (Color online). The minimum conductivity own 
versus the charged impurity distance d to the TI surface for 
nimp = 4.0 x 10 13 cm -2 . The solid, dashed and dot-dashed 
lines correspond to the TFDT-EMT, quasi-TFDT and the 
self-consistent methods, respectively, (a) The effective back- 
ground dielectric constant k = 5. (b) The effective back- 
ground dielectric constant re = 50. 
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C. Transport at finite temperature 

In this section, using the 2-fluid model, we present our 
results for the conductivity in TIs at finite temperature. 
If we neglect the contribution of activation processes in 
TIs the dominant contribution to the temperature depen- 
dence of a is due to electron-phonon scattering processes. 
Fig. [14] shows the longitudinal acoustic phonon limited 
resistivity of Bi 2 Se3 surface as a function of temperature 
on a log- log plot. To calculate the Bi2Se3 surface re- 
sistivity due to phonon scattering^ 7 , we use D = 30 eV 
for the deformation potential coupling constant, p m ~ 
7.68 x 10~ 7 g/cm -2 for the two dimensional mass den- 
sity, i.e. 1 quintuple layer mass density of Bi2Se3, and 
vi = 2900 m/s is the velocity of the longitudinal acoustic 
phonon mode~. The inset of Fig. [14] shows the log- 
arithmic derivatives of the temperature-dependent resis- 
tivity. Fig. [14] clearly demonstrates two different regimes 
depending on whether the phonon system is degener- 
ate or non-degenerate, and the low- to high-temperature 
crossover is characterized by the Bloch-Gruneiscn (BG) 
temperature Tbg = 2&fVz / fcg2£ r 37 i 64 . The resistivity in- 
creases with T as p ~ T A at low temperatures and p ~ T 
at high temperatures, which agrees with the results ob- 
tained for TI films by using an isotropic elastic contin- 
uum approach^ 7 -. We note that electron-phonon scat- 
tering is an important scattering mechanism for finite- 
temperature transport in TIs, e.g., Tbg ~ 30 K in TIs 
compared to Tbg ~ 100 K in graphene because the lat- 
ter has much larger phonon velocity. In graphene, in 
contrast to 2D surface TI transport, phonon effects are 
extremely weak^. 




FIG. 14. (Color online). Double-logarithmic scale of resis- 
tivity (Bi2Se3) only due to phonon scattering as a function 
of temperature for several densities. The deformation po- 
tential coupling constant D — 30 eV, p m — 7.68 x 10~ 7 
g/cm~ 2 , and the velocity of longitudinal acoustic phonon 
mode vi = 2900 m/s. The inset shows the logarithmic deriva- 
tives dlogp/dlogT versus temperature. 
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FIG. 15. (Color online). Calculated conductivity as a func- 
tion of temperature for different carrier density including only 
charged impurity scattering with d = 0.1 nm, « = 50 and 
nimp = 4 x 10 13 cm" 2 , (a) The charged impurity distance 
d = 0.1 nm. The corresponding TFDT n rms ^ 2.27 x 10 12 
cm -2 and the potential fluctuation s ~ 220 meV (Eq. ll4[) . (b) 
The charged impurity distance d = 1 nm. The corresponding 
TFDT n rms ~ 1.0 x 10 12 cm" 2 and the potential fluctuation 
s ~ 150 meV (Eq.[[4]). 



The results shown in Fig. [13] do not include the con- 
tribution to the resistivity due to quenched disorder. To 
include both the effect of quenched disorder and electron- 
phonon scattering we use the quasi-TFDT+2-fhiid. As 
discussed in Sec. HH this approach allows us to take into 
account several finite temperature effects: the temper- 
ature dependence of the screening of the quenched dis- 
order, electron-phonon scattering processes, broadening 
of the Fermi surface, and temperature induced activated 
processes. The temperature activated processes cause 
a to increase with T and therefore induce an insulat- 
ing behavior (i.e. conductivity increasing with increas- 
ing temperature) for cr(T), whereas the electron-phonon 
scattering processes induce a metallic behavior (i.e. con- 
ductivity decreasing with increasing temperature). The 
change with T of the screening of the disorder potential 
also induces a metallic behavior for cr(T), however given 
the large value of K in typical TIs this effect is quite 
weak, contrary to the case of graphene 10 or 2D semicon- 
ductor systems^. Figs. [1~5] (a), (b) show the dependence 
of a on T obtained by neglecting the effect of electron- 
phonon scattering processes. We see that in this case 
the temperature dependence of a is almost completely 
determined by thermally activated processes that induce 
a monotonically increase of a with T . 

Fig. HH shows the scaling of the conductivity with re- 
spect to doping at different temperatures including both 
the effects of quenched disorder and electron-phonon 
scattering. At large densities the main effect of the fi- 
nite temperature is to suppress a due to the presence of 
electron-phonon scattering processes. However, at low 
densities the effect of electron-phonon scattering pro- 
cesses competes with thermal activation processes and 
can give rise to a non-monotonic dependence of a with re- 
spect to T. This is shown in Figs. 1171 (a), (b) where <r(T) 
for different values of n is plotted. From Figs. [17] (a), (b) 
we see that for n < n rms at low temperatures the thermal 
activation processes dominate and induce an insulating 
behavior for <r(T). The crossover temperature from in- 
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sulating to metallic behavior depends on n, n n np, « and 
d. In general the larger the strength of the spatial fluc- 
tuations of the carrier density the stronger is the effect 
of thermal activation processes and therefore the larger 
is the low temperature range for which the transport ex- 
hibits insulating behavior. This is shown clearly by the 
scaling of <r(T) at the CNP for different values of ni mp , k, 
Figs. US] (a), (b), respectively. From Figs. [T51 (a) and (b) 
we see that a change of n{ mp and k that increases n rms 
extends the range of temperatures over which cr m j n (T) 
exhibit an insulating behavior. 

The temperature dependence of a for different values of 
d is shown in Fig. [18] (c). It appears to contradict the gen- 
eral rule that a parameter change that increases n rms will 
increase the range of temperatures over which a(T) ex- 
hibit an insulating behavior. This is due to the combina- 
tion of two effects: (i) The fact that at large d (d> 1 nm) 
n rms is very low and so the resistivity pdis due to the 
quenched disorder at low carrier densities is much higher 
than the resistivity p p h due to electron-phonon scatter- 
ing, so that a tot (T)/<j tot (T = 5K) « cr dls (T) / a dls (T = 
5K); (ii) The decrease at large d of the metallic screening 
effects. The scaling of a with respect to T and d is there- 
fore very interesting because it reveals the temperature 
dependence of the screening and could therefore be used 
to indirectly identify the nature of the disorder potential 
and the screening properties of the surfaces of 3DTIs. 
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FIG. 17. (Color online). Calculated conductivity as a func- 
tion of temperature for different carrier density including 
both charged impurity and electron-phonon scattering, with 
d = 0.1 nm, k = 50 and n imp = 4 x 10 13 cm -2 , (a) The 
charged impurity distance d 
TFDT n rms ~ 2.27 x 10 12 cm 
s ~ 220 meV (Eq. 1 14|) . (b) The charged impurity distance 
d = 1 nm. The corresponding TFDT n rms ^ 1.0 x 10 12 cm -2 
and the potential fluctuation s ~ 150 meV (Eq. I14[). 
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FIG. 16. (Color online). Calculated conductivity as a func- 
tion of net carrier density for various temperatures with d = 
0.1 nm and k = 50. (a) The impurity density n,i mv = 4 x 10 13 
cm" 2 . The corresponding TFDT n rms 2.27 x 10 12 cm" 2 
and the potential fluctuation s ~ 220 meV (Ea.lT4"]l; (b) The 
impurity density ni mp = 11 x 10 13 cm" 2 . The corresponding 



TFDT n rms ~ 4.21 x 10 1 
s ~ 300 meV (Eq. [HJ 
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Before concluding this section on the temperature- 
dependent surface conductivity of 3DTIs, we point out 
that one of the important qualitative findings of our work 
is the nonmonotonic temperature dependence of the 2D 
surface conductivity as apparent in Figs. ll6lll8l and as ex- 
pected from the competing mechanisms of phonon scat- 
tering and disorder induced density inhomogeneity. In 
particular, phonons induce higher-temperature metallic 
temperature dependence and the density inhomogene- 
ity induces insulating temperature dependence through 
thermal activation, and at some disorder-dependent (and 
also doping-dependent) characteristic temperature the 



FIG. 18. (Color online). Minimum conductivity as a func- 
tion of temperature with (a) for various charged impurity 
density n,i mp . As ni mp increases, the corresponding s = 
220, 300, 360, 390 meV; (b) for various dielectric constant 
k. As k increases, the corresponding s = 840, 500, 220, 110 
meV; (c) for various charged impurity distance d. As d in- 
creases, the corresponding s = 220, 180, 150, 130 meV. The 
corresponding potential fluctuation s is obtained via Eq. [14] 
using n rms calculated within the TFDT (shown in Figs. [3] [S] 
and©. 



transport behavior changes from being insulating-like 
to metallic-like. We emphasize that this temperature- 
induced crossover behavior has nothing to do with any 
localization phenomenon (and in fact, Anderson local- 
ization effects are completely absent in TI surface trans- 
port since all back scattering is suppressed) , and it arises 
entirely from a competition between inhomogeneity and 
phonons. The metallic behavior moves to higher (lower) 
temperature as disorder increases (decreases). This non- 
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monotonicity has been observed in experiments at lower 
2D carrier densities 3 ^ where the inhomogeneity effects 
are important. 

IV. CONCLUSIONS 

We have presented a detailed theoretical study of the 
2D transport properties of the surfaces of 3DTIs. There 
is compelling evidence that in current transport exper- 
iments on 3DTIs charged impurities are the dominant 
source of quenched disorder. For this reason in our 
study we have considered in detail the case in which the 
quenched disorder potential is the one created by ran- 
dom charged impurities close to the surface of the TI. 
However, the theoretical framework that we have devel- 
oped and presented in this work is very general. As in 
graphene, the presence of charged impurities induces the 
formation of strong carrier density inhomogeneities. In 
particular, close to the charge neutrality point the car- 
rier density landscape breaks-up in electron-hole puddles. 
The strong carrier density inhomogeneities make the the- 
oretical description of the electronic transport challeng- 
ing for two reasons: (i) inability to use standard theoret- 
ical methods that assume a homogeneous density land- 
scape; (ii) the importance at finite temperature of ther- 
mally activated processes. Our work presents a theoreti- 
cal description of transport on the surface of 3DTIs that 
overcomes these difficulties and takes into account also 
the effect of electron-phonon scattering processes which 
is an important resistive mechanism at higher tempera- 
tures. 

To characterize the disorder-induced inhomogeneities 
we use the Thomas-Fermi-Dirac-Theory. We also present 
a comparison of the results obtained using the TFDT 
with the one obtained using the self-consistent approx- 
imation and the quasi-TFDT approach and show that 
the three approaches give results that are qualitatively 
similar but that differ quantitatively. 

To study the electronic transport in the presence of 
strong carrier density inhomogeneities, starting from the 
TFDT results, we use the effective medium theory and a 
2-fluid model. The TFDT+EMT approach is well justi- 
fied and is expected to provide the most accurate results. 
However, the generalization of the TFDT+EMT method 
to finite temperature is impractical due to the contribu- 
tion to transport of thermally activated processes, con- 
tribution that can be dominant at low temperatures due 
to the strong carrier density inhomogeneities. At finite 
temperature the 2-fluid approach is very valuable be- 
cause it allows to take into account all the finite tem- 
perature effects, such as temperature dependent screen- 
ing and electron-phonon scattering, including thermally 
activated processes. The 2-fluid model relies on the use 
of an effective Gaussian distribution for the probability 
distribution of the screened potential P(V). The param- 
eters that define P(V) can be chosen in such a way to 
maximize the agreement of the results for T = obtained 



using the 2-fluid model and the ones obtained using the 
TFDT+EMT approach. These parameters are then used 
to obtain the transport properties at finite temperature 
using the 2-fluid model. 

In current 3DTIs the dielectric constant (k ~ 50) is 
much larger than in graphene where n ~ 1 — 4. In ad- 
dition, in 3DTIs the acoustic phonon velocity is smaller 
than in graphene. These facts make the contribution 
of electron-phonon scattering processes to the resistiv- 
ity much more important in the surface of 3DTIs than in 
graphene. For the surfaces of 3DTIs the effect of electron- 
phonon scattering events becomes important already for 
T as low as 10 K, whereas in graphene it becomes relevant 
only for T > 200 K. As a consequence for the surfaces of 
3DTIs we find that electron-phonon scattering is much 
more important to determine the dependence of the con- 
ductivity on T. The large value of k in 3DTIs also implies 
that for the surface of 3DTIs, contrary to graphene, the 
temperature dependence of the screening does not play 
an important role. The temperature dependence of the 
conductivity on the surface of 3DTIs is therefore mostly 
determined by electron-phonon scattering processes and 
thermal activations processes. These two types of pro- 
cesses have opposite effect on a and at low temperature 
and low doping compete giving rise to a nonmonotonic 
dependence of a with respect to T. The nonmonotonic 
temperature dependence of the 2D surface conductivity 
is one of the important new qualitative results of our 
theory. We have presented detailed results for a(T) that 
clearly show the competition of the different processes 
that affect a for T ^ 0. 

The theoretical approach developed here, being able 
to include all the main effects that determine the trans- 
port properties of the surfaces of 3DTIs, allowed us to 
present results that can be directly and quantitatively 
compared to the experimental ones. The good agree- 
ment between our theoretical results and the recent ex- 
perimental measurements ^ 15 ' 34 , suggest that the theo- 
retical method presented is very effective to characterize 
the transport properties of the surfaces of 3DTIs espe- 
cially due to its ability to take into account the effects due 
to the disorder-induced carrier density inhomogeneities. 
Future improvement of the theory could include a more 
accurate surface band structure and the effects of the 
bulk bands, but we do not expect these details to affect 
our qualitative results at low surface doping densities be- 
cause our theory includes the most important resistive 
processes contributing to the surface transport in 3DTIs. 
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